Beyond the Neutron Drip-Line: The Unbound Oxygen Isotopes 25 and 2e O 
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The very neutron-rich oxygen isotopes 25 O and 26 O are investigated experimentally and theoret- 
ically. In this first R3B-LAND experiment, the unbound states are populated at GSI via proton- 
knockout reactions from 26 F and 27 F at relativistic energies around 450 MeV/nucleon. From the 
kinematically complete measurement of the decay into 24 O plus one or two neutrons, the 25 O ground- 
state energy and lifetime are determined, and upper limits for the 26 O ground state are extracted. 
In addition, the results provide evidence for an excited state in 26 O at around 4 MeV. The ex- 
perimental findings are compared to theoretical shell-model calculations based on chiral two- and 
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three-nucleon (3N) forces, including for the first time residual 3N forces, which are shown to be 
amplified as valence neutrons are added. 



PACS numbers: 21.10.-k, 25.60.t, 27.30. +t, 29.30. Hs 

Understanding the properties of nuclei with extreme 
neutron-to-proton ratios presents a major challenge for 
rare-isotope beam experiments and nuclear theory. Nu- 
clei located at and beyond the neutron dripline play a 
crucial role in this endeavor. Experimentally, the neu- 
tron dripline has been established up to oxygen pQ with 
24 being the last bound isotope, while it extends consid- 
erably further in fluorine [2] , providing first experimental 
evidence for the particle stability of 31 F and particle in- 
stability of 28 0. Recently, it has been shown that the 
anomalous behavior in the oxygen isotopes is due to the 
impact of 3N forces, which provide repulsive contribu- 
tions to the interactions of valence neutrons [3] , connect- 
ing the frontier of neutron-rich nuclei to the theoretical 
developments of nuclear forces. 

Another striking feature in oxygen is the doubly magic 
nature of 22 and 24 [TJ 3HS] m strong contrast to the 
lighter elements, where the dripline is marked by nuclei 
exhibiting a loosely-bound halo structure. The neutron- 
rich oxygen isotopes also provide interesting insights, 
when viewed coming from stable nuclei. As protons are 
removed along isotones, the attractive contribution from 
the proton-neutron tensor force decreases, thus, opening 
up the N = 16 neutron shell for oxygen [T, which is 
below the N — 20 magic number present in stable nuclei. 

How the structure proceeds beyond 24 towards N = 
20 is thus of central interest. The present experimen- 
tal limits are 25 - 26 0, for which a resonance energy and 
width have been reported for 25 O [4] and a resonance 
position for 26 O [8]. Those experiments had a limited 
acceptance for neutrons, resulting in a narrow energy 
window to which the experiment was sensitive [S]. In 
this first R3B-LAND experiment, we investigate both un- 
bound isotopes 25 > 26 with an essentially flat efficiency 
up to 8 MeV relative energy. It has been speculated that 
the unbound isotopes 26 and 28 might have a rather 
long lifetime, which would constitute the first example 
of neutron radioactivity [TO]. Our present result estab- 
lishes a firm upper limit for the lifetime of the 26 O ground 
state. We then combine the experimental investigation 
with theoretical calculations based on chiral two-nucleon 
(NN) and 3N forces, where we focus on the increasing 
contribution from residual three-neutron forces as neu- 
trons are added. 

The experiment has been carried out at the GSI 
Helmholtzzcntrum in Darmstadt using the R3B-LAND 
reaction setup. Beams of light neutron-rich nuclei 
have been produced by fragmentation reactions of a 
490 MeV/nucleon 40 Ar primary beam in a 4 g/cm 2 Be 
target. Ions with a magnetic rigidity of 9.88(±1%) Tm 
corresponding to an A/Z ratio of around 3 have been se- 
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FIG. 1. (color online) Relative-energy spectrum of the 
24 0+n system measured in a proton-knockout reaction from 
26 F. The dotted curve shows a Breit-Wigner fit to the data. 
The grey histogram indicates the experimental response for 
a white spectrum. The inset displays x 2 as a function of 
position and width, and indicates the minimum and the la- 
contour. 

lected by the fragment separator FRS and transported 
to the experimental area. Energy-loss and time-of-flight 
measurements served for identification of the incoming 
ions on an event-by-event basis. The beam cocktail con- 
tained 26 > 27 F ions, which were selected to populate the 
unbound states in 25 . 26 0, by one-proton knockout reac- 
tions. Different targets have been used, while most of 
the statistics stems from measurements with the CH2 
(922 mg/cm 2 ) target. The target was surrounded by the 
47T detector Crystal Ball consisting of 160 Nal crystals for 
detecting photons and light particles emitted at large an- 
gles. Position and energy loss of the beam and fragments 
behind the target were detected by two silicon strip de- 
tectors before deflection in the large-gap dipole magnet 
ALADIN. Two further position measurements after the 
magnet with scintillating fiber detectors served for track- 
ing of the ions through the dipole field. Together with 
time-of-flight and energy-loss measurements, this yields 
the magnetic rigidity and atomic number and thus the 
mass of the fragments. Neutrons from the decay of un- 
bound states were detected at a distance of around 12 m 
downstream the target by the LAND neutron detector 
with an efficiency of 92% for one neutron events and cov- 
ering an angular range of ±80 mrad around the beam 
axis. A similar experimental setup and analysis scheme 
is described in Ref. [11] in more detail. 

From the measurements of the momenta of outgo- 
ing fragments and neutrons, the two- and three-body 
relative-energy spectra are reconstructed for one- and 
two- neutron events. Figure [l] shows the 24 0+n relative- 
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FIG. 2. (color online) Width and lifetime as a function of 
resonance energy for 26 O (upper panel) and 25 O (lower panel). 
The curves show theoretical expectations for different /-values 
of the neutron(s) from Ref. [10]. For 26 0, the upper limits 
for the resonance energy and lifetime are given by horizontal 
and vertical grey lines. The allowed region defined by these 
limits is represented by the shaded area. For 25 0, the present 
experimental results are given by the horizontal and vertical 
grey lines, with errors indicated by the shaded area. 



energy spectrum after proton removal from 26 F. A promi- 
nent peak structure is visible at around 700 keV corre- 
sponding to the ground-state resonance of 25 O. The po- 
sition and width of the resonance have been extracted 
by fitting a Breit-Wigner distribution with an energy- 
dependent width and a channel radius of 4 fm. For 
the x 2 minimization procedure, Pearson's chi-square 
method [T2], using errors of the parent distribution ac- 
cording to a Poisson probability distribution, has been 
used, as the usual method with errors estimated from the 
number of counts gives inaccurate results in case of low 
statistics. The lcr-contour in the 2-dimensional parame- 
ter space is displayed in the inset of Fig. [T] The extracted 
position and width of 737±H keV and 72±lf keV, re- 
spectively, are in good agreement with the previously re- 
ported values [3] within experimental errors, see Table [I] 
The relatively large experimental error on the width is 
due to the instrumental energy resolution which domi- 
nates the apparent width. The measured width is also 
in agreement with the expectation of a d-state single- 
particle resonance with a spectroscopic factor of unity, as 
shown in the lower panel of Fig.[2]by the expected widths 
and lifetimes as a function of resonance position for dif- 
ferent neutron angular momentum I (based on Ref. |10) ) . 
The experimental relative-energy spectrum for 26 is 
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FIG. 3. (color online) One- and two-neutron channel mea- 
sured in coincidence with 24 O after one-proton removal from 
27 F. The histograms with error bars represent the measured 
24 0+2n three-body relative-energy spectrum (top panel) and 
the 24 0+n+X relative-energy spectrum (bottom panel). The 
grey histogram indicates the experimental response for a 
white spectrum in the 2n channel. The dotted curves dis- 
play the most probable fit to the data corresponding to two 
states in 26 at 25i 2 j keV and around 4 MeV (see text). 



shown in Fig. [3j The upper panel displays events, where 
24 has been detected in coincidence with two neutrons. 
Two groups of events are observed: below 1 MeV and 
around 4 MeV. The shape of the experimental response, 
indicated by the grey histogram, is rather flat over the 
energy region shown, but exhibits a steep fall-off for small 
energies below 500 keV. For such small relative energies, 
the neutrons are not well separated in space and time 
when hitting the detector and can thus hardly be distin- 
guished from In events. The effect can be seen quanti- 
tatively in the two-dimensional response matrices shown 
in Fig. [4] The energy reconstructed from the simulation 
is plotted versus the generated one in the lower panel, 
showing a band along the diagonal with a width reflecting 
the instrumental resolution. For low relative energies, the 
curve bends away from the diagonal and, in addition, the 
events are being reconstructed as In events with a large 
probability. This can be seen in the upper panel, which 
shows the reconstructed 24 0+n relative energy spectrum 
for the events falsely identified as In events. The simula- 
tion is based on measured real In events from deuteron 
breakup reactions. The 2n events are generated by over- 
laying the hit patterns from these measured In events 
and are then analyzed in the same way as the experi- 
mental data. The generated response matrices thus do 
not rely on a simulation of the neutron reactions in the 
detection process. 

The effect discussed above is clearly visible in the 26 O 
data of Fig. [3] In the In channel (lower panel), the ac- 
cumulation of events at very low energy in the first bin 
is striking. This is not present in the In events mea- 
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FIG. 4. (color online) Simulated instrumental response for 
detection of a 2n decay of 26 O with relative energy E re i. The 
lower panel shows the extracted energy for two detected neu- 
trons in coincidence with 24 O, while in the upper part, the 
24 0+n relative energy spectrum is shown for events, where 
only one neutron has been detected due to inefficiency of the 
detector or the multi-neutron recognition algorithm. 



TABLE I. Compilation of experimental and theoretical re- 
sults obtained in this work compared to previously published 
results. All energies and lifetimes are given in keV and ns, re- 
spectively. Energies are measured from the gs energy of 24 0. 



Er r t Ref. 

25 0(gs) 737+1 72t£f 9T±f -lCT 12 this work 

770l 2 g 172±jg Bl 

733 NN+3N + residual 3N this work 

1301/1303 USDA/B pj] 

1002 - - [T5] 

2e O(gs) ^40/1200 ^1.2 10" 1( |] sJ5.7[] this work 

150+^p : m 

18 NN+3N + residual 3N this work 

501/356 USDA/B [14] 

21 0^02 - [15] 

26 0(exc) 4225l 22 6 this work 

a 68%/95% c.l. 

b from lifetime estimate, see text and Fig. 2 
c 95% c.l. 



sured in the proton knockout from 6 F, see Fig.^ The 
events in the first bin of the O+n spectrum are the 
characteristic fingerprint of a very low-lying state in 26 O. 
We have performed a simultaneous statistical analysis of 
In and 2n coincidences, starting with the hypothesis of 
a low-lying state in 26 0. The dotted histogram displays 
the most probable result yielding a position for the 26 O 
ground state of 25+ 2 g keV. Again, the \ 2 minimization 
using Poisson distributed errors of the test function has 
been used. 

An alternative method |13j . which is independent of 
the binning of the experimental data, is to use for each 
event the measured relative energy in conjunction with 
the inverse response matrix to calculate the probability 
distribution of the true energy (Bayesian approach with 
uniform prior). The resulting Bayesian interval goes from 
to 40/120 keV for 68%/95% confidence limit (c.l.), in- 
cluding the In events at very low relative energy. From 
the energy measurement alone, we cannot exclude a very 
small value close to zero for position and width poten- 
tially leading to a rather long-lived 26 ground state, 
which would constitute the first case of neutron radioac- 
tivity as speculated in Ref. 10 : As can be seen from 
Fig. [2J the lifetime for a <i 2 -state could reach seconds for 
a resonance position well below 1 keV. Such a long life- 
time, however, can be excluded from the fragment mea- 
surement. The distance of the target to the middle of the 
dipole magnet measures 256 cm corresponding to a flight 
time of 11.8 ns. If 26 decayed after that time, a frag- 
ment mass larger than 24 would be reconstructed by the 
tracking procedure. Also, no neutron coincidences should 
be observed in that case, since the fragment is bent by 
7 degree after passing half of the dipole field. From the 



analysis of the energy spectra discussed above, the num- 
ber of produced 26 can be estimated. By triggering 
on proton knockout reactions using the Crystal Ball and 
vetoing on the neutron detector, the fragment-mass spec- 
trum has been reconstructed for the 27 F(p, 2p) A reac- 
tion. While we still observe events corresponding to 24 0, 
only one event appears above the 24 cut at A = 24.8, 
which could potentially stem from a late decay in the 
dipole field. From this we obtain an upper limit for the 
lifetime of 5.7 ns with a 95% confidence limit. The upper 
limit for the position of the state and for its lifetime are 
shown in Fig. [2] leaving the shaded area as the allowed 
region. There is only a small overlap with the calcu- 
lated values for a cP-state, suggesting a more complex 
26 O ground-state configuration. 

It would therefore be very interesting to determine the 
lifetime experimentally more precisely in order to gain 
insights to the structure of 26 O. Such a measurement 
would be rather straightforward using a similar method 
as described here, but placing the target directly in front 
of the dipole. The decay curve of 26 O would translate 
into a fragment-mass spectrum distribution depending 
on the decay position. In addition, the neutrons will be 
detected not only at zero degree, but also in the bend- 
ing direction, again directly reflecting the decay curve. 
With the intensities available, e.g., at the PJBF facility 
at PJKEN, a precise value could be obtained from such 
an experiment with the SAMURAI detector. 

We compare the ground-state energies of 25 > 26 to the- 
oretical shell-model calculations based on chiral effective 
field theory interactions combined with renormalization- 
group methods to evolve nuclear forces to lower resolu- 
tion [TB]. Our results are based on chiral NN and 3N 
forces, where the single-particle energies and two-body 
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FIG. 5. (color online) Comparison of the experimental 
25 O and 26 O energies with theoretical shell-model calcula- 
tions based on chiral NN and 3N forces (NN+3N) and includ- 
ing residual 3N forces. The impact of residual 3N forces is 
highlighted by the arrows. 



interactions of valence neutrons on top of a 16 core are 
calculated following Refs. [171II8], without adjustments. 
Here, we focus on the relative contribution from resid- 
ual three-valence-neutron forces, which become more im- 
portant with neutron number along isotopic chains [T§j . 
Figure [5] shows the predicted ground-state energies. The 
relative contribution from residual 3N forces, which we 
included perturbatively, increases from 0.1 to 0.4 MeV for 
24-26Q Because the ground states of 25 ' 26 are very nar- 
row, thus quasi-bound, they can be treated fairly well in a 
bound-state approximation. The remarkable agreement 
with experiment should, however, be taken with care, 
as the continuum needs to be included. The expected 
contribution is 200 keV to both 25 ' 26 [2U], so relative 
energy differences will be smaller. We have also explored 
uncertainties in the calculation of the single-particle en- 
ergies and valence interactions, which would increase 
the ground-state energy relative to 24 for 25 > 26 by 
0.2 — 0.3 MeV. In Table|TJ we also compare with the phe- 
nomenological USDA/B interactions [TJj, which predict 
too high energies for 25 > 26 0. Better agreement is found in 
Table [i] for the continuum shell model [15] , see also the 
promising continuum developments of Refs. [20"1 121j. 

The theoretical calculations based on chiral NN and 
3N forces predict a first excited 2 + state in 26 O at 
1.6 MeV relative to the 26 ground state. It is found at 
1.9/2.1 MeV for USDA/B interactions and at 1.8 MeV 
for Ref. [TS] • Experimentally, the events at 4 MeV in the 
three-body energy spectrum (top panel of Fig. [3]) provide 
evidence for an excited state in 26 O, with most proba- 
ble energy 4225t 22 £ keV. A likely candidate would be a 
proton-hole state populated after knockout from the 27 F 
Opi/2 shell (rather than from the 0d 5 / 2 valence orbital). 
To investigate proton (and neutron) cross-shell excita- 



tions, we have carried out (0 + l)hu! calculations in the 
spsdpf space using the WBP interaction [22 , for which 
the first and second 2 + energies are at 2.3 and 5.4 MeV, 
repectively. The first state with a proton excitation com- 
ponent from 0pi/2 to 0d 5 /2 is a 3" state at higher en- 
ergy of 5.4 MeV, whose proton contributon is also mixed 
with neutron excitations and considerably weaker than 
for the corresponding 3~ state at 5.0 MeV in 18,20 O. Note 
that 26 O is bound by 1.0 MeV for the WBP interaction, 
so that 1 MeV uncertainties are possible. The lowest 
negative parity states predicted by the WBP interaction 
are a quartet of 3~,2~, 1 — ,0 _ neutron excitations from 
0g? 3 /2 to lp 3 / 2 at 3.7, 4.1, 4.5, 4.9 MeV (with centroid at 
4.1 MeV). Because of strong proton- neutron interactions, 
the proton-knockout reaction mechanism could also be 
more involved, generating admixed neutron excitations. 

In this first R3B-LAND experiment, we have investi- 
gated the ground-state energy, width, and lifetime of the 
unbound oxygen isotopes 25 O and 26 O. Our results are 
in very good agreement with theoretical shell-model cal- 
culations based on chiral NN and 3N forces, where the 
ground-state energy of these extremely neutron-rich iso- 
topes becomes increasingly sensitive to 3N forces among 
the valence neutrons. The 26 ground state is unbound 
by less than 120 keV, and our measurement for the life- 
time ^5.7 ns (both at 95% c.l.) excludes 26 being a 
first candidate for neutron radioactivity. We have pro- 
posed an experiment with existing detectors to determine 
the 26 O lifetime accurately, which would yield important 
information on the ground-state nature of 26 0. We also 
obtained first evidence for an excited state of 26 O located 
at around 4 MeV. Different possibilities for the nature of 
this state exist, making it an exciting case for future cal- 
culations and experiments with higher resolution. 
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